A low level of extragalactic background light 
as revealed by /-rays from blazars 
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The diffuse extragalactic background light consists of the sum of the starlight emitted by galaxies 
through the history of the Universe, and it could also have an important contribution from the first 
stars, which may have formed before galaxy formation began. Direct measurements are difficult and 
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not yet conclusive, owing to the large uncertainties caused by the bright foreground emission associ- 
ated with zodiacal light 1 . An alternative approach 2-5 is to study the absorption features imprinted 
on the 7-ray spectra of distant extragalactic objects by interactions of those photons with the back- 
ground light photons 6 . Here we report the discovery of 7-ray emission from the blazars 7 H 2356— 309 
and 1ES 1101—232, at redshifts z=0.165 and z=0.186, respectively. Their unexpectedly hard spectra 
provide an upper limit on the background light at optical/near-infrared wavelengths that appears to 
be very close to the lower limit given by the integrated light of resolved galaxies 8 . The background 
flux at these wavelengths accordingly seems to be strongly dominated by the direct starlight from 
galaxies, thus excluding a large contribution from other sources - in particular from the first stars 
formed 9 . This result also indicates that intergalactic space is more transparent to 7-rays than previ- 
ously thought. 

The observations were carried out with the High Energy Stereoscopic System 10 (H.E.S.S. ), a system 
of four imaging atmospheric Cherenkov telescopes operating at energies E > 0.1 TeV. These two blazars 
are at present the most distant sources for which spectra have been measured at these energies (Tab. 1). 

Intergalactic absorption is caused by the process of photon-photon collision and pair production. The 
original spectrum emitted by the source (which we call "intrinsic") is modified such that the observed flux 
Fohs{E) = F- mt (E) ■ e~ T ( E \ where the optical depth r(E) depends on the Spectral Energy Distribution 
(SED) of the Extragalactic Background Light (EBL) (Fig. 1). Details are provided in the Supplementary 
Notes and Figures. For any reasonable range of fluxes at ultraviolet (UV) and optical/near-infrared wave- 
lengths (O-NIR), t(E) - and thus absorption - is larger at 1 TeV with respect to 0.2 TeV. This difference 
makes the observed spectrum steeper (that is, r o b s > Ti nt , for a power-law model dN/dE oc E~ r ) The 
spectral change Ar=r o b s — Tint scales linearly with the EBL normalization, and becomes more pronounced 
at larger redshifts. Thus more distant objects provide a more sensitive diagnostic tool. 

In general, if the intrinsic spectrum were sufficiently well known, t(E) - and thus the EBL SED - 
could be effectively measured by comparing intrinsic with observed spectra. Blazars, however, are char- 
acterized by a wide range of possible spectra, and the present understanding of their radiation processes is 
not yet complete enough to reliably predict their intrinsic 7-ray spectra. But for these two sources, with 
O-NIR fluxes at the level of the "direct" estimates, the intrinsic spectra needed to reproduce the H.E.S.S. 
data become extremely hard (that is, they have small values of T), at odds with the currently known blazar 
physics and phenomenology. This can be avoided by reducing the energy dependence of the optical depth, 
assuming either very low O-NIR fluxes (reducing r at 1 TeV) or very high UV-0 fluxes (increasing r at 
0.2 TeV). The latter case, however, would require unreasonably high UV fluxes (as discussed later). We can 
then derive an upper limit on the EBL by requiring the intrinsic spectrum to be compatible with the present 
knowledge of blazars. 

To determine such a limit, a plausible shape for the EBL SED is assumed. As a reference shape in 
the 0.1-10 pLva. range we adopted the phenomenological curve used in refs 6 ' 11 , which is designed to be in 
general agreement with the EBL spectrum expected from galaxy emission 12 ' 13 . This curve, labelled P1.0 in 
Figs 1 and 2, was originally normalized to match the "direct" estimates at 2.2 and 3.5 jtxm (refs 1 ' 14 ' 15 ). Here 
we leave its normalization as a free parameter, scaling P 1.0 by different factors P (labelled accordingly, the 
curve scaled by 0.45 x is "P0.45") down to the lower limit obtained by the resolved galaxy counts 8 (~P0.4). 
To reproduce the excess around 1.5 fim claimed from the Infrared Telescope in Space (IRTS) data 16 (also 
argued by ref 17 ), an additional ad hoc component was considered, labelled "Enir". This feature is not 
expected from standard galaxy evolution models, and could be the spectral signature of radiation produced 
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in the early universe, for example by the first stars formed (metal-free massive stars, called 'Population III'; 
see e.g. refs 9 ' 18 ). 

The intrinsic source spectra (Fig. 2) have been reconstructed directly from the observed ones using 
the assumed EBL, without a priori assumptions on the blazar spectrum. EBL evolution effects (e.g. due to 
galaxy evolution) were not included: these effects, negligible at low redshifts, become important as redshift 
increases, but for the range considered here (z = 0.165 — 0.186) their impact is still limited to a factor of 

< 10% (AP < 0.1; see Supplementary Information). The reconstructed spectra are generally compatible 
with a power-law (dN/dE oc E~ r ), but the EBL densities PI.O+Enir and P1.0 both yield extremely hard 
photon indices (r; nt < 0, see Fig. 2 and Supplementary Fig. SI-2), implying a pile-up or line-like feature 
in the blazars' SED around 1-3 TeV. We obtain the same result by considering the NIR excess added to the 
galaxy counts limits (e.g. P0.4+Enir, see Supplementary Fig. ST3). This is because a lower EBL flux only 
in the UV-0 band decreases the absolute values of r but increases the contrast between 0.2 and 1 TeV. 

Such hard spectra have never been seen in the closest, less absorbed TeV blazars such as Mkn 421 and 
Mkn 501 19 ~ 21 (z ~ 0.03, r int « 1.5 - 2.8 using the same EBL SEDs), and are difficult to explain within 
the present standard leptonic or hadronic scenarios 6 for blazar emission. In shock acceleration models, 
the hardest index obtained for the accelerated particles is s = 1.5 (see e.g. ref 22 ). In the case of protons 
interacting with ambient plasma, the resulting 7-ray spectrum has the same slope, r; nt = 1.5. In the 
case of electrons, the spectrum of the 7-rays emitted through inverse Compton scattering is expected to 
be steeper than T^t = 1.5 under most circumstances. Only if radiative cooling is not effective and the 
blazar Compton emission is wholly within the Thomson limit - unlikely at such high energies - one finds 
Tint = (s+ 1)/2 = 1.25. We thus assume in the following discussion that the true average intrinsic spectrum 
was not harder than r; nt = 1.5, although later we also address the possibility of harder spectra. 

To be compatible with r; nt > 1.5, the EBL flux has to be scaled down to P0.45 to explain both 
objects' data, with 1ES 1101—232 providing the most stringent constraints thanks to the better statistics at 
high energies and the larger redshift. With a fixed EBL shape, there is a direct link between the normalization 
P and Tint. The one-sigma statistical and systematic uncertainties on the H.E.S.S. spectral measurement can 
then be translated to an equivalent uncertainty on P, AP ~ ±0.15 (see Supplementary Information). 

This limit (P0.45) is robust with respect to a different EBL spectral shape, as long as it maintains an 
overall maximum around 1-2 fim. Below 1 /im, lower fluxes than our template tend to harden the intrinsic 
spectrum more, while even a flat slope from 1.4 fim down to 0.1 /jm would soften it by only AT ~0.1. 
Above 2 /im, the slope cannot be much flatter than our template - a flatter slope could explain in part the 
"direct" measurements at 3.5 //m (Supplementary Fig. SI-4) - because this would again imply a new, very 
hard component (r; nt < 0) in the intrinsic spectrum, rising at a few TeV (Supplementary Fig. SI-5). In this 
respect, this H.E.S.S. dataset gives the same indication as the HEGRA data 11 on H 1426+428 (z=0.129), 
which show a flattening feature above 1 TeV naturally provided by a starlight EBL between 3 and 10 /jm 
(SED oc A" 1 ). 

Therefore, the conservative and self-consistent assumptions of both not unusual blazar spectra (Tint > 
1.5) and a galaxy-like EBL spectrum allow the EBL flux around 1-2 fim to be constrained at the level of 

< (14 ± 4)nW m~ 2 sr^fhat is, < 0.55 ± 0.15 xPl.0). This corresponds to P(0.45+0.1) to allow for 
galaxy evolution effects. Coupled with the lower limits derived from galaxy counts given by the Hubble 
Space Telescope 8 (~9.0-9.7±f;g nW m 2 sr 1 ), the H.E.S.S. spectra lead to conclude that more than 2/3 
of the EBL in the O-NIR band is resolved into single sources. This result is completely independent of 
any "direct" measurement of the EBL. Remarkably, it is in severe conflict with the claims of high EBL flux 



3 



at NIR wavelengths ' and, to a lesser extent, with the reported detections at 2.2 and 3.5 /j,m * . The 
H.E.S.S. upper limits agree instead with the most recent theoretical calculations 23 of the EBL, as well as 
with recent theoretical arguments 24 ' 25 against high EBL fluxes due to Pop III stars. 

This result is also rather insensitive to small changes in the assumed intrinsic slope. A different value, 
if proved more likely by future results, will shift the limit accordingly, but only strong spectral differences 
would qualitatively change our conclusion: even a value of T- mt = 1.0 would loosen the 0.55 x limit only to 
< 0.7xP1.0. 

Alternative scenarios which could reconcile the measured spectra with high O-NIR fluxes formally 
exist, and would represent a major discovery in their own right, but we consider them very unlikely, given 
their exotic implications. Higher UV fluxes would make the intrinsic spectra softer, but the huge values 
required (> 300 nW m~ 2 sr _1 , see Fig. 1, for example) are not supported by other measurements 1 ' 26 , and 
could not be accomodated within any reasonable cosmological model. 

A more viable alternative is that such hard spectra are a real, new feature of the TeV blazar emission. 
Possible mechanisms have already been envisaged 6 . For example, the inverse Compton scattering of mono- 
energetic electrons (Eq, such as cold plasma with very large bulk motion Lorentz factor), interacting in deep 
Klein-Nishina regime with a narrow-band photon field (such as a Planck-type distribution) may lead to very 
flat 7-ray spectra with a sharp pile-up at e 7 ~ Eq, reproducing spectra like the ones in Fig. 2. However, such 
features should become directly visible in the observed spectra of the closer, less absorbed objects of the 
same type, like the well-studied Mkn 421 and Mkn 501 (if r; nt was ks 0, they should show r o b s < 0.5). This 
is in contrast with observations 19 21 , unless we assume a dependence of the source parameters on redshift 
such that the corresponding features always disappear due to EBL absorption. It is difficult to justify such a 
fine-tuning on a relatively small redshift range, although more objects and observations are needed to settle 
this issue definitively, given the still-limited sample. 

Other possibilities include the non-cosmological origin of blazars' redshifts and the violation of 
Lorentz invariance (see e.g. ref 27 ). However, both scenarios imply dramatic revisions of modern physics 
and astrophysics, which we do not consider to be justified by these data alone. 

A low EBL level, in agreement with the expectations from standard galaxy evolution models, is the 
simplest and most likely explanation of the H.E.S.S. data. 
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1ES1101-232 


0.186 


43 


~12 


0.16-3.3 


2.88 ±0.17 


4.44 ± 0.74 • 10" 


-13 


0.62/11 


H 2356-309 


0.165 


40 


~10 


0.16-1.0 


3.06 ±0.21 


3.08 ±0.75 • 10" 


-13 


0.66/6 



Table 1 : Main parameters of the H.E.S.S. observations. These observations were performed in 
June-December 2004 for H 2356-309, and March-June 2004 & 2005 for 1ES1101-232. The 
table gives the total exposure after selection for good quality data, significance of the detected 7- 
ray signal, energy range used for the spectral fits and the result of a single power-law fit (dN/dE = 
N (E/TeV)~ T ° hs ). The spectra have been calculated applying the technique described in ref 10 . 
Errors are \o statistical. The systematic uncertainty on the flux and photon index are estimated 
to be ~15% and ~0.1, respectively. Details of these observations will be published elsewhere; 
here we have focused on the cosmological implications of the measured spectra. Compared to 
the previous observations of TeV blazars, these H.E.S.S. spectra provide significantly stronger 
constraints on the EBL because of the combination of a hard spectrum and relatively high redshift 
(see Supplementary Information). 
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Figure 1 Spectral Energy Distribution of the EBL in the wavelength band most affecting these 
H.E.S.S. data (0.1-10 /mi). The EBL data are from a review compilation 1 (errors ^a), unless other- 
wise stated. Open symbols correspond to the integrated light from galaxy counts, and thus must 
be considered lower limits for the EBL: in the UV-0 range, from Hubble data (green, red 8 ); in 
the NIR, from Spitzer (blue 28 ) and ISO data. Note that these data are also lower limits for the 
total emission from galaxies, because of various observational and selection effects in the detec- 
tion and counting of faint galaxies. The possible missed light in the the UV-0 band has been 
estimated 29 to be < (2-3) nW m~ 2 sr _1 . The upper limits (purple) are 2a estimates 1 . Direct mea- 
surements are shown as filled symbols: IRTS data from the NIR Spectrometer 16 (blue), and data 
from COBE/DIRBE (green 15 , magenta 17 and red triangles). Red squares correspond to tentative 
detections in the optical 26 with corrections according to ref 30 . The curves show the EBL shapes 
used to reconstruct the intrinsic spectra. P1 .0 gives 26, 23 and 14 nW m 2 sr x at 1 .25, 2.2 and 
3.5 fim, respectively. The thick line shows the range most effectively constrained by the H.E.S.S. 
data. In the long-dashed regions, higher fluxes than P0.55 would not be in conflict, as long as the 
flux in the 1-3 /mi range is within or around the limit. The short dashed line shows the additional 
UV component needed by P1 .0 to soften the intrinsic spectrum down to r = 1.5 (see Supplemen- 
tary Fig. SI-5; E N i R would require even higher fluxes). This example is the most energetically 
economic solution, limited to the narrow range ~ 0.2 - 0.4 /xm to have the maximum effect on the 
7-ray spectrum with the minimum UV flux and the minimum impact on the overall attenuation. 
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Figure 2 The H.E.S.S. spectra of 1ES 1101-232, corrected for absorption with three different 
EBL SEDs, as labelled in Fig. 1. Red: observed data. Blue: absorption-corrected data. The 
data points are at the average photon energy in each bin, also used to calculate the optical 
depth for reconstruction. For the calculation, a flat A-dominated cosmology was adopted, with 
H = 70 km/s/Mpc, Vt m = 0.3, fi A = 0.7. Error bars are 1a, statistical errors only. Between 
1 .3 and 3.3 TeV, the overall detection is 4a. The lines show the best-fit power laws to the recon- 
structed spectrum (dN/dE = N (E/TeV)~ Tint ), and the corresponding shapes after absorption. 
The Xrec/d-o.f. (calculated by integrating the absorbed power-law model over each observed data 
bin) are (from left to right): 1.20/11, 0.54/11, 0.47/11. We note that possible spectral variability 
does not weaken our conclusions because it would imply states with even harder spectra than the 
average one (by definition). Note also that the X-ray spectrum (which in blazars usually samples 
the synchrotron emission of TeV electrons) measured during simultaneous observations in June 
2004 and March 2005 does not show such hard slopes, but is similar to the historical values (Aha- 
ronian et al. 2006, in preparation). For H 2356-309, the same EBL SEDs yield r int = -0.6, 0.7 
and 2.0, respectively (Suppl. Fig. SI-2). The NIR excess onto the galaxy counts limits (P0.4+E N i R ) 
yields r int « -0.7 and -0.4 for the two objects, respectively (see Suppl. Fig. SI-3). 
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Supplementary Information 



In this paper we have derived an upper limit on the Extragalactic Background Light (EBL) in the optical to 
near infrared band (O-NIR). This is done by examining the effect of absorption, caused by interaction with 
EBL in this band, on the 7-ray spectra of two newly-detected distant blazars: 1ES 1101—232 (z=0.186) and 
H 2356—309 (z=0.165). Here we provide additional details on the manner in which EBL absorption modifies 
the source spectrum, and on the impact on the derived limit of the uncertainties in the 7-ray measurements 
and of galaxy evolution effects. We also explain why previous TeV blazar detections did not provide a 
comparably strong limit, and present supplementary figures to better illustrate the cases described in the 
paper. At the end, we provide a table with the H.E.S.S. spectral data for the two sources. 

1 Supplementary Notes 

EBL absorption 

Absorption modifies the intrinsic source spectrum according to F ^, S (E) = Fi nt (E) ■ e~ r( - E \ The optical 
depth t(E) is obtained by convolving the EBL photon density distribution n(e) with the pair production 
cross-section a(E, e), and integrating over the 7-ray photon path from the source to the observer: 



(where dl(z) = cdt(z) is the proper line element, a(E, e) is already integrated over angles for an isotropic 
field of background photons, and the photon density and energies depend on redshift). The energy depen- 
dence of t(E) strongly depends on the spectrum of the background photons (see e.g. ref[19]). For example, 
if n(e) oc e^ 1 (or equivalently the EBL SED vl v oc A -1 ), the optical depth appears constant with energy, 
while for n(e) oc e~ 2 (i.e., vl v = constant), the optical depth t{E) oc E. Note that in the first case, the 
observed spectrum simply reproduces the spectrum emitted by the source, attenuated by a constant factor. 
The energy dependence of the optical depth translates to strong modifications of the source spectrum even 
with relatively small differences in r, due to the exponential dependence. For example, given an intrinsic 
power-law spectrum dN/dE oc E~ Vint , a difference of only one unit in r between 0.2 and 1 TeV increases 
the spectral index by AT = r o b s — Ti nt = 0.62. Figure SI-1 shows the attenuation factor as a function of 
energy in four relevant cases (P1.0, P0.4, P0.4+E NIR , P1.0+UV). 

For a 7-ray photon of a given energy E, because of the threshold effect (Ee > (m e c 2 ) 2 ), only EBL pho- 
tons of energy above threshold (i.e., A < A t hr) contribute to the absorption. Also, because of the relative 
narrowness of the 77 — > e + e~ cross-section, for broad-band EBL spectra more than half of the interactions 
occur with a narrow interval of target photons AA ~ (1 ± 1/2)A* (ref[6,19]) around the cross-section peak, 
where A* « 1.4 {E 1 / 1 TeV) /im (this relation is illustrated on the upper axis in Fig. SI-3 and SI-4, as a 
guideline). Thus, the optical depth at ~ 0.2 TeV is controlled only by the EBL fluxes at UV-0 wavelengths 
(below 1 micron), since NIR photons are below the energy threshold, while the optical depth at ~ 1 TeV is 
controlled mainly by the NIR flux. The UV-0 band has a more limited impact on ~ 1 TeV photons since the 
interactions occur far from the cross-section peak. An increase or decrease of the UV-0 flux alone changes 
the optical depth at 0.2 TeV much more than at 1 TeV. 

Therefore, it is possible to reduce the energy dependence of the optical depth (and thus to avoid the hard 
reconstructed spectra) by assuming either very low O-NIR fluxes (reducing r at TeV energies) or very high 
UV-0 fluxes (increasing r at sub-TeV energies; e.g., see Fig. SI-1, right panel). 




(1) 
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Comparison with earlier observations 

These two new objects provide significantly stronger constraints on the EBL than any previous TeV blazar 
because of a favourable combination of hard observed spectrum and relatively high redshift. For sources at 
smaller redshifts Ar is smaller, while at the same redshift a softer spectrum can be misinterpreted as a hard 
spectrum which has been more heavily absorbed. Both cases allow the EBL to have a wider range of possible 
fluxes without yielding implausible values in the reconstructed blazar spectrum. This is the reason why the 
TeV spectra of the two previous most distant objects, 1ES 1426+428 [11] (z=0.129, r obs ~ 3.5 ± 0.4) 
and PKS 2155-304 [10] (z=0.116, r obs ~ 3.3 - 3.4 ± 0.06), did not provide such strong constraints as 
these two new sources, even when measured with much higher signal-to-noise ratio (PKS 2155— 304). The 
intrinsic spectra of H 1426+428 and PKS 2155-304 start to become very hard (T < 1) for EBL fluxes 
around P1.0+Enir and above. 

Statistical and systematics errors 

The determination of an EBL limit from gamma-ray spectra needs to take into account the statistical and 
systematic errors of the spectral determination, which amount for these data - at the la level - to AT ~ 0.2 
and 0.1, respectively. In addition, the absolute energy scale of Cherenkov telescopes is only determined 
to ~ 15%. The EBL limit can be different, therefore, because for example the true blazar spectrum can 
be softer than the H.E.S.S. value, and the photon energies lower. Using the same Tint = 1.5 limit on the 
intrinsic spectrum, these errors can be translated into an equivalent uncertainty for the scaling factor P: 
namely, by finding the values for which T^t ± Ari ntjS tat ± Ar^sys = 1-5, and shifting the bin energies 
±15%. For the 1ES 1101—232 spectrum, which provides the most stringent constraints, the errors in F 
translate to AP ~ ±0.1 (the relation is linear, and is approximately AP ~ 0.34AT), and the energy scale 
error into AP ~ ±0.05. All three contributions (±lcr statistical, ±0.1 systematic, ±15% energy scale 
uncertainty) combined yield AP ~ ±0.15. 

Galaxy evolution effects 

In the calculation of the optical depth, we have not included the effects of galaxy evolution between z=0 and 
z = 0.165 to 0.186. That is, we assume a constant (rather than decreasing with redshift) EBL co-moving 
energy density. This corresponds to the assumption that all the background photons seen today were already 
in place at the source redshift. To estimate the impact of this approximation in the redshift range considered 
here, we used the model by Primack [12] in 2001, interpolating a set of EBL SEDs as a function of redshift 
(kindly provided by the authors). Including galaxy evolution, the intrinsic spectra become AT < 0.2 softer 
than without evolution, thus loosening the limit on the scaling factor P by less than 0.1. Note, however, 
that this effect can still be safely ignored when testing scenarios where the EBL is dominated by additional 
radiation produced in the early universe (for which our assumption fully holds). 
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Figure SI-1 Attenuation factor e~ r for the redshift of 1 ES 1 101-232 (z=0.186), assuming different EBL 
SEDs. Left panel: P1.0, P0.4, P0.4 with E Nm added on (Fig. SI-3). Right panel: P1.0 and P1.0+UV. 
These curves would correspond to the observed TeV spectrum for a source with a flat power-law emission 
spectrum (e.g., for a source spectral energy distribution E 2 dN/dE with r int = 2). Above 0.2 TeV, the 
increase of r with E steepens the spectrum, up to 1-2 TeV. Left panel: the thick lines (corresponding 
to P0.4, P1.0) show how an increase of the EBL normalization steepens the 0.2-1 TeV spectrum. The 
dashed line (P0.4+E Nm ) shows how the addition of a NIR excess around ~ l(j,m dramatically steepens the 
spectrum, due to the increase of r for 0.8-1 TeV photons, while photons at 0.1-0.3 TeV remain unaffected. 
Above 2 TeV, the absorption curve corresponding to P0.4+E N i R even hardens with energy since P0.4+E Nm 
is characterized by a slope steeper than A -1 between 2 and 6 fim. The right panel shows how an increase 
in the UV flux reduces the energy dependence of the optical depth, at the price of a higher attenuation 
overall. 
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H 2356-309 




Figure SI-2 The H.E.S.S. time-averaged spectrum of H 2356-309 (green), and the absorption-corrected 
data (blue) with the same EBL SEDs as for 1 ES 1 1 01 -232 (Fig. 1 and 2). The data points and the 99% 
upper limit are at the average photon energy for each bin, also used to calculate the optical depth for recon- 
struction. Error bars are \a statistical errors. The lines show the best-fit power laws to the reconstructed 
spectrum (dN/dE = N (E/TeV)~ Fint ), and the corresponding shapes including absorption. The x? od /d.o.f. 
are (from left to right): 1 .94/6, 0.86/6, 0.59/6. 
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Figure SI-3 Effect of the NIR excess onto the galaxy counts limits. Left: EBL SED with the same data 
as Fig. 1. The full line shows the EBL shape with the NIR excess added on to P0.4, while the right 
panel shows the corresponding absorption-corrected 7-ray spectrum for 1 ES 1 1 01 -232. Also, for this EBL 
SED, the intrinsic spectrum is extremely hard (r int = -0.7 between 0.1 and 2 TeV), with a very sharp rise 
between 0.3 and 1.5 TeV where the intrinsic slope becomes as hard as r int = -1.2. This is due to the 
strong difference in optical depths at these two energies, because the gamma-ray photons below 0.3 TeV 
are insensitive to the NIR excess > 1/xm (since they are below the energy threshold for the 77 — > e + e~ 
process), whereas ~ lTeV photons are fully affected since the interaction takes place near the peak of 
the 77 — > e+e~ cross-section. The exclusion of the the IRTS data points between 1 .4 and 2.2 jum (dashed 
line) is not enough to avoid very hard spectra (r int « 0.2): significantly lower NIR fluxes are needed, as 
described in the main text. We note that ref[25] argues that the shape of the NIR excess between 1 and 
4/im is almost identical to the spectrum of the zodiacal light, and could therefore be of local instead of 
extragalactic origin. 



13 



0.22 0.73 2.2 



EjJeV) 



100 




Pl.O 



P0.55 



0.1 1 10 

X[/u.m] 

Figure SI-4 Effect of modifications to the EBL SED in the UV and IR bands, for the same data as 
Fig. SI-2. The short dashed line shows the additional UV component needed by P1 .0 to soften the intrinsic 
spectrum down to r int = 1.5, as shown in Fig. SI-2. This is the most energetically-economic solution, limited 
to the narrow range ~ 0.2-0.4 to have the maximum effect on the 7-ray spectrum with the minimum UV 
flux and overall attenuation. The inclusion of E Nm would require even higher UV fluxes, in the case of both 
P1.0+E N i R and P0.4+E Nm . Above 2 /mi, a slope much flatter than our template (e.g. +IR, dot-dashed line) 
is not supported by either the H.E.S.S. (see Fig. SI-5) or HEGRA results [11], which instead suggest an 
EBL spectrum around oc A -1 . If this is the case, it would constrain the EBL to be near the recent Spitzer 
source counts [28], meaning that we are close to resolving the diffuse background also up to ~ 8/xm. 
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Figure SI-5 The H.E.S.S. spectrum of 1 ES 1 1 01 -232, corrected for absorption with the EBL SED shown 
in Fig. SI-4. The P1.0+UV panel shows the spectral softening given by the additional UV component, at 
the price of increased energetic requirements: the intrinsic TeV flux would be 3 orders of magnitude higher 
than that observed, and ~ 1.5-2 orders of magnitude above the blazar synchrotron peak, as indicated by 
simultaneous X-ray data (Aharonian et al. 2006, in preparation). The P0.55+IR panel shows the hardening 
given by an EBL spectrum above 2/im which is flatter than our template (Fig. SI-4, dot-dashed line). 
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Table S 1-1 : H.E.S.S. spectral data. Tha table gives the differential flux $ in different energy bins, for 
both objects. E avg is the average photon energy energy in each bin. The upper limit for H 2356-309 
is given for the overall energy range [1.041 - 3.292] TeV, with a confidence level of 99%. 



Energy Bins 




1ES 1101-232 








H 2356-309 




Elow Efrigh 


E av g 


$± A$ 




Eavg 


$± A$ 




[TeV] 


[TeV] 


[cm- 2 s- 1 TeV- 1 




[TeV] 


[cm^s^TeV- 1 




0.165 - 0.208 


0.184 


5.38 ±1.88 


x 10" 


il 


0.184 


4.58 ±1.34 x 10' 


-ll 


0.208 - 0.262 


0.232 


2.93 ±0.78 


x 10" 


-li 


0.232 


3.08 ±0.66 x 10' 


-li 


0.262 - 0.329 


0.292 


1.67 ±0.39 


x 10' 


-li 


0.292 


1.37 ±0.33 x 10' 


-li 


0.329 - 0.414 


0.367 


9.00 ±2.08 


x 10' 


-12 


0.367 


5.40 ±1.57 x 10' 


-12 


0.414 - 0.522 


0.462 


4.57 ±1.04 


x 10' 


-12 


0.462 


4.60 ±0.94 x 10' 


-12 


0.522 - 0.657 


0.582 


1.41 ±0.54 


x 10' 


-12 


0.582 


1.54 ±0.52 x 10' 


-12 


0.657 - 0.827 


0.733 


1.17 ±0.37 


x 10' 


-12 


0.732 


5.94 ±2.96 x 10' 


-13 


0.827 - 1.041 


0.922 


4.22 ±2.34 


x 10' 


-13 


0.921 


3.64 ±1.89 x 10' 


-13 


1.041 - 1.311 


1.161 


1.59 ±1.13 


x 10' 


-13 








1.311 - 1.650 


1.462 


1.83 ±1.02 


x 10' 


-13 


1.572 


<5.03 x 10- 14 




1.650 - 2.077 


1.840 


1.17 ±0.56 


x 10' 


-13 








2.077 - 2.615 


2.316 


8.49 ±3.87 


x 10' 


-14 








2.615 - 3.292 


2.916 


4.73 ±3.03 


x 10- 


-14 
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